SPOT: Sliced Partial Optimal Transport

International audienceOptimal transport research has surged in the last decade with wide applications in computer graphics. In most cases, however, it has focused on the special case of the so-called ``balanced'' optimal transport problem, that is, the problem of optimally matching positive measures of equal total mass. While this approach is suitable for handling probability distributions as their total mass is always equal to one, it precludes other applications manipulating disparate measures.Our paper proposes a fast approach to the optimal transport of constant distributions supported on point sets of different cardinality via one-dimensional slices. This leads to one-dimensional partial assignment problems akin to alignment problems encountered in genomics or text comparison. Contrary to one-dimensional balanced optimal transport that leads to a trivial linear-time algorithm, such partial optimal transport, even in 1-d, has not seen any closed-form solution nor very efficient algorithms to date.We provide the first efficient 1-d partial optimal transport solver. Along with a quasilinear time problem decomposition algorithm, it solves 1-d assignment problems consisting of up to millions of Dirac distributions within fractions of a second in parallel.We handle higher dimensional problems via a slicing approach, and further extend the popular iterative closest point algorithm using optimal transport -- an algorithm we call Fast Iterative Sliced Transport. We illustrate our method on computer graphics applications such a color transfer and point cloud registration

Efficient Distance Transformation for Path-based Metrics

In many applications, separable algorithms have demonstrated their efficiency to perform high performance volumetric processing of shape, such as distance transformation or medial axis extraction. In the literature, several authors have discussed about conditions on the metric to be considered in a separable approach. In this article, we present generic separable algorithms to efficiently compute Voronoi maps and distance transformations for a large class of metrics. Focusing on path-based norms (chamfer masks, neighborhood sequences...), we propose efficient algorithms to compute such volumetric transformation in dimension $n$. We describe a new $O(n\cdot N^n\cdot\log{N}\cdot(n+\log f))$ algorithm for shapes in a $N^n$ domain for chamfer norms with a rational ball of $f$ facets (compared to $O(f^{\lfloor\frac{n}{2}\rfloor}\cdot N^n)$ with previous approaches). Last we further investigate an even more elaborate algorithm with the same worst-case complexity, but reaching a complexity of $O(n\cdot N^n\cdot\log{f}\cdot(n+\log f))$ experimentally, under assumption of regularity distribution of the mask vectors

Optimal Separable Algorithms to Compute the Reverse Euclidean Distance Transformation and Discrete Medial Axis in Arbitrary Dimension

In binary images, the distance transformation (DT) and the geometrical skeleton extraction are classic tools for shape analysis. In this paper, we present time optimal algorithms to solve the reverse Euclidean distance transformation and the reversible medial axis extraction problems for $d$-dimensional images. We also present a $d$-dimensional medial axis filtering process that allows us to control the quality of the reconstructed shape

Variance Analysis for Monte Carlo Integration: A Representation-Theoretic Perspective

In this report, we revisit the work of Pilleboue et al. [2015], providing a representation-theoretic derivation of the closed-form expression for the expected value and variance in homogeneous Monte Carlo integration. We show that the results obtained for the variance estimation of Monte Carlo integration on the torus, the sphere, and Euclidean space can be formulated as specific instances of a more general theory. We review the related representation theory and show how it can be used to derive a closed-form solution

Dynamic Reconstruction of Complex Planar Objects on Irregular Isothetic Grids

International audienceThe vectorization of discrete regular images has been widely developed in many image processing and synthesis applications, where images are considered as a regular static data. Regardless of final application, we have proposed in [14] a reconstruction algorithm of planar graphical elements on irregular isothetic grids. In this paper, we present a dynamic version of this algorithm to control the reconstruction. Indeed, we handle local refinements to update efficiently our complete shape representation. We also illustrate an application of our contribution for interactive approximation of implicit curves by lines, controlling the topology of the reconstruction

Discrete bisector function and Euclidean skeleton in 2D and 3D

International audienceWe propose a new definition and an exact algorithm for the discrete bisector function, which is an important tool for analyzing and filtering Euclidean skeletons. We also introduce a new thinning algorithm which produces homotopic discrete Euclidean skeletons. These algorithms, which are valid both in 2D and 3D, are integrated in a skeletonization method which is based on exact transformations, allows the filtering of skeletons, and is computationally efficient

Integral based Curvature Estimators in Digital Geometry

International audienceIn many geometry processing applications, the estimation of differential geometric quantities such as curvature or normal vector field is an essential step. When designing such estimators, we have to pay attention to both its theoretical properties and practical effectiveness. In this paper, we investigate a new class of estimators on digital shape boundaries based on Integral Invariants. More precisely, we provide proofs of multigrid convergence of curvature estimators which are easy to implement on digital data. Furthermore, we discuss about some algorithmic optimisations and detail a complete experimental evaluation

Implementation of Integral based Digital Curvature Estimators in DGtal

In many geometry processing applications, differential geometric quantities estimation such as curvature or normal vector field is an essential step. In [1], we have defined curvature estimators on digital shape boundaries based on Integral Invariants. In this paper, we focus on implementation details of these estimators

Scale-space Feature Extraction on Digital Surfaces

International audienceA classical problem in many computer graphics applications consists in extracting significant zones or points on an object surface,like loci of tangent discontinuity (edges), maxima or minima of curvatures, inflection points, etc. These places have specific localgeometrical properties and often called generically features. An important problem is related to the scale, or range of scales,for which a feature is relevant. We propose a new robust method to detect features on digital data (surface of objects in Z^3 ),which exploits asymptotic properties of recent digital curvature estimators. In [1, 2], authors have proposed curvature estimators(mean, principal and Gaussian) on 2D and 3D digitized shapes and have demonstrated their multigrid convergence (for C^3 -smoothsurfaces). Since such approaches integrate local information within a ball around points of interest, the radius is a crucial parameter.In this article, we consider the radius as a scale-space parameter. By analyzing the behavior of such curvature estimators as the ballradius tends to zero, we propose a tool to efficiently characterize and extract several relevant features (edges, smooth and flat parts)on digital surfaces

Integral based Curvature Estimators in Digital Geometry

International audienceIn many geometry processing applications, the estimation of differential geometric quantities such as curvature or normal vector field is an essential step. When designing such estimators, we have to pay attention to both its theoretical properties and practical effectiveness. In this paper, we investigate a new class of estimators on digital shape boundaries based on Integral Invariants. More precisely, we provide proofs of multigrid convergence of curvature estimators which are easy to implement on digital data. Furthermore, we discuss about some algorithmic optimisations and detail a complete experimental evaluation